o 



Quantum phase transitions in a new exactly solvable quantum spin biaxial model with 

multiple spin interactions 

A.A. Zvyagini'2 

^ Max- Planck- Institut fiir Physik komplexer Systeme, 

Nothnitzer Str., 38, D-01187, Dresden, Germany 

^B.I. Verkin Institute for Low Temperature Physics and Engineering of the National 

Academy of Sciences of Ukraine, Lenin Ave., J^l, Kharkov, 61103, Ukraine 

(Dated: July 14, 2010) 

The new integrable quantum spin model is proposed. The model has a biaxial magnetic anisotropy 
of alternating coupling between spins together with multiple spin interactions. Our model gives 
the possibility to exactly find thermodynamic characteristics of the considered spin chain. The 
ground state of the model can reveal spontaneous values of the total magnetic and antiferromagnetic 
moments, caused by multiple spin couplings. Also, in the ground state, depending on the strength 
of multiple spin couplings, our model manifests several quantum critical points, some of which are 
governed by the external magnetic field. 
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Integrable models of quantum physics of magnetism 
are, unfortunately, rare. Such models, however, are very 
important for theorists, because they permit to compare 
the results of standard for theoretical physics of real sys- 
tems perturbative approaches with exact ones. For ex- 
ample, the seminal Ising model served as a basis for many 
powerful methods of modern physics, like the scaling, 
renormalization group, etc. Quantum integrable mod- 
els in one space dimension (ID) are developed relatively 
more, comparing to the 2D or 3D ones (in fact there exist 
only few examples of quantum integrable models in 2D 
and 3D), due to the relative simplicity of their study. On 
the other hand, according to the Mermin- Wagner theo- 
rem, any nonzero temperature in ID (and 2D) destroys 
long-range magnetic ordering. That theorem reveals, in 
fact, the enhancement of quantum and thermal fluctua- 
tions in low-dimensional systems, due to peculiarities in 
their densities of states. 

On the other hand, the interest in quantum spin sys- 
tems, where spin-spin interactions along one or two space 
directions are much stronger than couplings along other 
directions, has considerably grown during last decade. 
Such interest to low-dimensional quantum spin systems 
is motivated, first of all, by the progress in the prepara- 
tion of real substances with well defined one-dimensional 
subsystems. On the other hand, modern technologies 
permit to compose artificial one-dimensional quantum 
systems, like quantum wires and rings, which properties 
are created to be similar to theoretically known mod- 
els. Devices, based on such especially prepared quantum 
ID spin systems are very promising in the modern nano- 
technology, in the development of spintronics, or even in 
the quantum computation. 

ID quantum spin systems often manifest quantum 
phase transitions, i.e. those, which take place in the 
ground state, and which are governed by other than the 
temperature parameters, like an external magnetic field, 
external or internal (caused by chemical substitutions) 
pressure, etc. 



In the past most of exactly solvable quantum spin mod- 
els were related to the class of models with only pair 
spin-spin interactions between nearest-neighbor spinsji 
Last years more attention of physicists was paid to the- 
oretical studies of quantum spin models with not only 
nearest-neighbor spin-spin interactions, but also with 
next-nearest neighbor ones, multiple spin exchange mod- 
els (e.g., a ring exchange), etc<^"— Such models some- 
times appear to be closer to the real situation in quasi- 
one-dimensional magnets, comparing to the ones with 
only nearest-neighbor couplings between only pairs of 
spins. For example, such additional interactions are often 
present in oxides of transition metals, where the direct 
exchange between magnetic ions is complimented by the 
superexchange between magnetic ions via nonmagnetic 
ones. Terms, involving the product of three and four 
spin operators or more, were only recently recognized 
to be important for the theoretical description of many 
physical systems, despite the fact that multiple spin ex- 
change models were introduced by Thouless already in 
1965.^^ For example, multiple spin exchanges were used 
for the description of the magnetic properties of solid 
•^He^iS. Later similar models were used to study some 
cupratesi^ and spin ladders J^ Often quantum spin mod- 
els with antiferromagnetic interactions and multiple spin 
interactions manifest the spin frustration, i.e. the low- 
est in energy state is highly degenerated^ Last but not 
the least, such models often reveal quantum phase tran- 
sitions. For instance, many transition metal compounds, 
like copper oxides, are believed to manifest features, char- 
acteristic for quantum phase transitions. However, it is 
known that for many quantum spin models the standard 
quasi-classical theoretical description, based on the quan- 
tization of small deviations of classical vectors of magne- 
tization of magnetic sublattices, yield incorrect results, 
especially in the vicinity of quantum critical points.^ This 
is why, quantum exactly solvable spin models with multi- 
ple spin exchange interactions, even being rather formal, 
and, sometimes, non-realistic, are of great importance: 



They provide the possibihty to check approximate theo- 
retical methods, used for the description of more reahstic 
physical models of quantum spin systems with spin frus- 
tration. 

In this paper we propose a new integrable model of 
quantum spins with nearest-neighbor interactions and 
multiple spin exchange. The aim of this work is to 
study a model, that, on the one hand, contains multiple 
spin interactions, which usually produce incommensurate 
magnetic structuresii Second, the proposed model con- 
sists alternating exchange interactions between nearest 
neighbors, which can be the reason for the spin gap for 
low-lying excitations;^ Finally, the model has the biax- 
ial magnetic anisotropy, which is believed to be the key 
property of transition metal compounds with strong spin- 
orbit coupling J^ The 2D counterpart of the model can be 
related to the plaquette model of pi -1- ip2 superconduct- 
ing arrays.— On the other hand, the model is relatively 
simple, because the Hamiltonian of the model can be ex- 
actly transformed to the one of the free fermion lattice 
gas, and, hence, most of thermodynamic characteristics 
can be calculated explicitly. 

The Hamiltonian of our exactly solvable quantum spin 
model with alternating nearest-neighbor couplings and 
three-spin interactions, which permits exact solution, has 
the form: 
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where S^'1'2 are the operators of the spin-^ projections 
of the spin in the n-th cell, which belongs to the sub- 
lattice 1 or 2, /ii,2 are effective magnetons of the sublat- 
tices, H is the external magnetic field, directed along z 
axis, Ji^2.x,y are the alternating exchange coupling con- 
stants between nearest neighbor spins in the cell and be- 
tween cells, and Ji,23 are alternating coupling constants 
for three-spin interactions. Notice that the model re- 
veals the biaxial magnetic anisotropy, i.e. the exchange 
interactions (in the spin subspace) along x, y, and z di- 
rections are different. In the case Ji^2x = Ji,2y, i-e. in 
the case of only uniaxial magnetic anisotropy, the model 
contains, as a special case, the model, studied in Ref. |la. 
On the other hand, the special case J13 = J23 = of the 
model is known for many years J^ Finally, in the absence 
of the magnetic field, H = 0, and three-spin couplings, 
Ji3 — J23 — the model can be related to the so-called 
ID quantum compass model (in the special case Jix — a, 
Jiy = 1 — a, J2X = J2y = 1)^ Terms, in the Hamilto- 
nian, which describe three-spin couplings, obviously vi- 
olate time-reversal and parity symmetries of the system 
separately, but the combination of both symmetries is 



preserved. 

After the Jordan- Wigner transformationii^ 
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where aJn i 2 and a„^i^2 are creation and destruction oper- 
ators, which satisfy fermionic anticommutation relations, 
and, after the Fourier transform 



an,i,2 = N ^^^ ^ afc,i,2 exp(z/cn) 



(3) 



and similar for a^ 121 where N is the number of cells, the 
Hamiltonian Eq. ([T]) gets the form 
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where J12 ~ (1/2)(Ji,2k ± Ji,2y)- With the help of a 
unitary transformation this Hamiltonian can be diago- 
nalized 



fc w— 1 V ^ 



where 



k j=i 



^U,2^Fk±dFi-Gk, 



(5) 



(6) 



and 



+ 2{JixJ2y + JlyJ2x) COS fc + ( jfg + jfg) COS^ fc 

+A{nl + nl)H'^ - A{^J^lJl^, + H2J23)H cosk 




3 3 




FIG. 1: Dispersion relations for the exactly solvable quantum 
compass model at zero magnetic field for the upper branch 
(grey) and lower branch (black) of eigenstates as functions of 
the parameter of three-spin interactions J23 for J13 — J23. 
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{Jix-hx - Jiyhy)'^ sin^ k 
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Using the standard particle-hole transformation one can 
get only positive eigenvalues of the Hamiltonian ([B]). One 
can check that in the case Ji,23 = the spectrum coin- 
cides with the one of Ref. [l^ while for Ji^2x = Ji.2y it 
reproduces the spectrum from Ref. [l^ The energies of 
eigenstates of the first branch are non-negative for all 
parameters of the model. The energies of eigenstates, 
belonging to the second branch can be equal to zero, de- 
pending on the values of the parameters of the model. 
Figs. [T][2] represent the zero field dispersion relations for 
both branches as functions of three-spin interactions for 
the homogeneous three-spin couplings and for the alter- 
nating three-spin couplings, respectively, for the quan- 
tum compass model with a = 0.4 

Despite some parameter-dependent features, the be- 
havior of our model is similar for the quantum compass 
case with very strong biaxial magnetic anisotropy and 
the case with small biaxial anisotropy. 

It is simple to obtain thermodynamic characteristics of 
our model at nonzero temperatures. The free energy of 
the quantum spin chain is equal to 
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FIG. 2: The same as in Fig. 1, but for J13 = 0. 



Obviously, the z-projection of the average magnetization 
of the system is 



AP 
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From this formula it is easy to show that M^ is zero 
for H = for any nonzero temperature, in accordance 
with the Mermin- Wagner theorem. The low temperature 
behavior of the magnetic susceptibility. 
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and the specific heat. 
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depend on the values of coupling constants J12, </i3,23, 
the effective magnetons, /ii,2j and the value of the ex- 
ternal magnetic field H , see below. One can check that 
there is no ordering and, therefore, none of thermody- 
namic characteristics of the considered system has pecu- 
liarities at any nonzero temperature. On the other hand, 
as it will be shown below, in the ground state sponta- 
neous magnetic ordering can take place. In the cases, 
where elementary excitations of the model are gapped, 
the low-temperature magnetic susceptibility and the spe- 
cific heat reveal exponential in T dependencies in the 
absence of spontaneous magnetization. If the model re- 
veals the spontaneous magnetic moment, the magnetic 
susceptibility at low temperatures is divergent. On the 



other hand, for gapless situation of fow-energy states of 
the model, the magnetic susceptibility is finite at low 
temperatures for the absence of spontaneous magnetic 
ordering at T = 0, while the specific heat is linear in 
T . At the critical lines of quantum phase transitions 
(see below) our model manifests either square root, or 
logarithmic in T and magnetic field behaviors of the spe- 
cific heat and the magnetic susceptibility. In the case, 
where interaction constants are very different from each 
other (or, to be more precise, when two branches of eigen- 
states are characterized by very different energy scales), 
the specific heat and the magnetic susceptibility can re- 
veal two-maxima temperature dependencies. 

The most important properties of the one-dimensional 
spin system are manifested in the ground state. The 
ground state energy of our model can be written as 



En = ^ 




(12) 



Then the z-projections of each total spin moment of a 
cell in the ground state can be written as: 
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The sum of the ^-projections of spin moments can be con- 
sidered as the ground state vector of magnetism, or mag- 
netization of the model, M^ = fiiSf + /J.2<5'|, while the 
difference describes the vector of antiferromagnetismfiS 
L^ = ^iS^ — 1^282 , or the staggered magnetization of the 
model. From these expressions one immediately sees that 
in the ground state the model can have nonzero spon- 
taneous magnetic and antiferromagnetic moments (i.e. 
magnetic ordering for iJ = 0), caused by nonzero three- 
spin interactions. We would like to turn attention that 
the signs of Ji,23 do matter. Namely, depending on their 
signs, the spontaneous magnetization of the model in the 
ground state can be positive or negative, with respect to 
the direction of the magnetic field. It is different from 
the behaviors of other known exactly solvable spin chain 
models. The reason for the onset of the spontaneous 
magnetic and antiferromagnetic moments in the ground 
state of our model is related to the violation of the time- 
reversal symmetry by three-spin coupling terms. 



It is interesting to notice that the equality Gfe = 
means that ek,2 = 0. As it is shown below, namely the 
condition G^ = determines the features in the behavior 
of all ground state characteristics of the spin chain. One 
can see, that Gk — either at sin fc = (i.e. for k — 0,Tr), 
or, for any k, if JixJ2x — JiyJ2y (it turns out that this 
condition does not depend on the magnetic field and on 
the values of three-spin couplings) . 

Let us consider first the case with JixJ2x 7^ JiyJ2y 
Notice that the limiting case of the quantum compass 
model belongs to the situation. The first branch of eigen- 
states is ever positive, but the second one can reach zero 
only for two values of the quasimomenta (k = 0, tt). Then 
it is simple to show that the ground state magnetisation 
is a continuous function of the external magnetic field, ex- 
cept of at i7 = for ^1 J23 ¥= ^2^13 and /ii J13 ^ M2./23, 
see Eqs. (|13m4p . For the latter the spontaneous magneti- 
zation appears, and, therefore, the ground state magnetic 
susceptibility is divergent there at H — 0. The magnetic 
susceptibility for nonzero values of H can have peculiar- 
ities, proportional to \n\H — Hc.i\ {i = 1,2,3,4), if the 
external magnetic field becomes equal to one of following 
four values 
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at which second order quantum phase transitions can 
take place, see Figs. [3]|31 

Such quantum critical points exist, naturally, only if 
those critical values of the field are real and non-negative. 
They are real if 

\^J■l^J■2{Jlx ± J2y){Jly ± J2x)\ < (^^1-^23 " M2>/l3)^ , (16) 

and the first two critical values are non-negative for 
/11M2 > 0, if 



(^^1^23 +M2J13) > 
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or the second two critical values are non-negative, if 
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For M1/X2 < the non-negativity conditions are reversed. 
If the reality conditions are not satisfied, no quantum 




FIG. 3: Critical values of the magnetic fields Hci,2 for the ex- 
actly solvable spin model as functions of parameters of three- 
spin interactions J13 and J23. We used ni = 1.01, IJ.2 ~ 0.99, 
Jix = 1, Jiy = 1.5, J2x = 2, and J2y = 0.9. 




FIG. 4: Critical values of the magnetic fields Hc3,4, for the ex- 
actly solvable spin model as functions of parameters of three- 
spin interactions J13 and J23. The set of parameters is the 
same as in Fig. 3. 



phase transitions, governed by the external magnetic 
field, take place in the system. If one of them is sat- 
isfied, and the other isn't, then only up to tvifo quantum 
phase transitions can happen. If the conditions Eqs. (|17p . 
or UHl), are not satisfied, then only one or two quantum 
phase transitions, governed by the field, take place. If 
one of the effective magnetons is zero (i.e. one of the 
ions, which form elementary cell, is non- magnetic), but 
three-spin interaction constants are not, then the quan- 



tum phase transitions take place at the values of the mag- 
netic field 



H. 



c5,6 



± 



J13J23 + {Jlx ± J2y){Jly ± J2x) 



2/^1,2-^; 



2,13 
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at which the magnetic susceptibility has logarithmic sin- 
gularities. Naturally, only positive values of iJc5,6 matter. 
Finally, if both of effective magnetons are zero, then, ob- 
viously, there are no quantum phase transitions, governed 
by the magnetic field. It turns out that at nonzero tem- 
peratures thermodynamic characteristics of the model re- 
veal logarithmic in T features at the critical values of the 
magnetic field. 

Consider now the situation, in which JixJ2x = JiyJ2y- 
The energies of the eigenstates in this case can be written 
as 

-1 2 
J23)cosfc 
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where 



A, 



Bk 



[Ji + J2 exp(-ifc)] 
[Jf - J2" exp(-ifc)] 



(21) 



It is obvious that Sk^i is positive for any parameters of 
the model. On the other hand, for the lower branch for 
some ranges of the quasimomentum k and external mag- 
netic field H, the first term under the square root sign in 
Eq. (|20l) can be smaller than the second one. It implies 
that eigenstates for lower branch can exist only for some 
ranges of fc, depending on the value of the external field 
H . The analysis of this situation is similar to the above 
(except the fact that one has to take into account nonzero 
Fermi seas, i.e. totally filled states with negative energies 
for some ranges of k depending on the value of the ex- 
ternal field; the critical value of k is determined from the 
condition e2,k^ = 0). One can see that there exist four 
critical values of the magnetic field, at which quantum 
phase transitions can take place, see Eqs. (ITS)) - ([T^ . The 
difference, comparing to the case with JixJ2x 7^ JiyJ2y, 
is in the more strong features of the magnetic suscep- 
tibility at critical fields ~ \/ yJ\H — Hc,i\ in the ground 
state, and, therefore, in square root peculiarities in T 
of thermodynamic characteristics of the model, like the 
magnetic susceptibility and specific heat, at critical val- 
ues of the magnetic field. 

Let us consider the homogeneous limiting case of our 
model Jix ^ J2x = Jx, Jiy = J2y = Jy, J13 = J; 
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Ml 
as 



M2 — M- In this case the Hamiltonian can be written 
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FIG. 5: The ground state dependencies of the magnetization 
as a function of the magnetic field for the homogeneous limit 
of the exactly solvable model for Jx — 1, Jy — 0.6 for J3 = —2. 



FIG. 6: The ground state dependencies of the magnetization 
as a function of the magnetic field for the homogeneous limit 
of the exactly solvable model for J^ = 1, Jy — 0.6 for J3 — 2. 



where 



fiH — -( J3 cos(2fc) + {Jx + Jy) cos(fc) 



".(^x Jy) 



Hk) 



(23) 



One can see that the energy (P5|) is non-negative. It can 



for 



be equal to zero only for J^ = Jy, or, if J^ 7^ Jy 
fc = 0, TT. In the later case there are two critical values of 
the magnetic field, at which quantum phase transitions 
can take place 



H 



{2tl)-'[J3±{Jx + Jy)] 



(24) 



Obviously, quantum phase transitions take place if values 
of the critical field are non-negative, i.e. they take place 

for ^ > 0, a 



J3 ± (J, + Jy)>0 



(25) 



Hence, J3 = ±( J^. + Jj,) are the conditions of the quan- 
tum phase transition, governed by the three-spin cou- 
pling. The ground state magnetic susceptibility has log- 
arithmic features ^ \n\h — H^i 2 1 ^^ critical values of the 
field. 

Figs. [S] and H] show the ground state behavior of the 
magnetization of our model for the homogeneous case as 
a function of the magnetic field. Fig. [5] presents the be- 
havior for J3 < —{Jx + Jy), and Fig. [HI demonstrates the 
magnetic field behavior for the region J3 > {Jx + Jy)- 
One can see that for both regions there is a spontaneous 
magnetization, but its sign (with respect to the direction 
of the field) depends on the sign of three-spin interac- 
tions. Also, there are two quantum phase transitions for 




FIG. 7: The ground state dependencies of the magnetization 
as a function of the magnetic field for the homogeneous limit 
of the exactly solvable model for J^ — 1, Jy = 0.6 for J3 = 0.2. 



J3 > {Jx + Jy), while for J3 < —{Jx + Jy) the ground 
state magnetization is a smooth function of H. 

Figs. [7] and [HI present the magnetic field behavior of the 
magnetization for —{Jx+Jy) < J3 < {Jx+Jy) for positive 
and negative values of J3 , respectively. One can see that 
in this region there is no spontaneous magnetization, and 
only one second order quantum phase transition takes 
place. 

On the other hand, if Jx — Jy, the eigenstates of the 
Hamiltonian can be negative for some ranges of k depend- 
ing on the value of the magnetic field. Negative energies 




H 



FIG. 8: The ground state dependencies of the magnetization 
as a function of the magnetic field for the homogeneous limit 



of the exactly solvable model for J^ 
-0.8. 



1, J„ = 0.6 for Ja 



imply the nonzero Fermi sea, where eigenstates with neg- 
ative energies are totally filled, and the ones with positive 
energies are empty. In that case quantum phase transi- 
tions yield square root singularities '^ 1/ J\H — Hj^i 2I 

of the magnetic susceptibility, cf. Refs. |^|^. 

It is important to point out that the quasiclassical de- 
scription of our model (when one replaces spin operators 
by classical vectors, and quantizing small deviations from 
the classical minimal energy state) does not reproduce ex- 



act results for the inhomogeneous (dimerized) situation. 
Namely, in the classical description of the model without 
biaxial anisotropy one of the branches of eigenstates is 
obviously gapless, unlike the exact result. 

In conclusion, motivated by recent experiments on 
quasi- ID quantum spin systems and recent theories of 
quantum compass model, we proposed the integrable 
model, in which exchange interactions between neighbor- 
ing spins is accompanied by the multiple spin exchange 
with the biaxial magnetic anisotropy. The model is sim- 
ple (due to the exact mapping to the problem of the lat- 
tice free fermion gas), and, therefore, permits to obtain 
exactly thermodynamic characteristics of the considered 
quantum spin chain. The most important behavior of 
the model is in the ground state. Our model manifests 
a ferrimagnetic-like ordering in the ground state. De- 
pending on the signs of the parameters of three-spin cou- 
plings, the spontaneous magnetic moment of the system 
in the ground state can be positive or negative (with re- 
spect to the direction of the magnetic field). The system 
can undergo several second order quantum phase tran- 
sitions, governed by the external magnetic field and the 
three-spin couplings strengths (the later can be caused 
by an external or internal pressure). Despite some artifi- 
cial structure of our model, we expect that more realistic 
quantum biaxial spin systems with multiple exchange in- 
teractions and the alternation of the exchange between 
nearest neighbor spins, will show similar to our simple 
model behavior, i.e. our exact solution has generic fea- 
tures for this class of quantum systems. 
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